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^ " ABSTRACT 

Context. Rotational splittings are currently measured for several main sequence stars and a large number of red giants with the space 
mission Kepler. This will provide stringent constraints on rotation profiles. 

Aims. Our aim is to obtain seismic constraints on the internal transport and surface loss of angular momentum of oscillating solar-like 
stars. To this end, we study the evolution of rotational splittings from the pre-main sequence to the red-giant branch for stochastically 
excited oscillation modes. 

Methods. We modified the evolutionary code CES AM2K to take rotationally induced transport in radiative zones into account. Linear 
rotational splittings were computed for a sequence of 1.3M Q models. Rotation profiles were derived from our evolutionary models 
and eigenfunctions from linear adiabatic oscillation calculations. 

Results. We find that transport by meridional circulation and shear turbulence yields far too high a core rotation rate for red-giant 
models compared with recent seismic observations. We discuss several uncertainties in the physical description of stars that could 
have an impact on the rotation profiles. For instance, we find that the Goldreich-Schubert-Fricke instability does not extract enough 
angular momentum from the core to account for the discrepancy. In contrast, an increase of the horizontal turbulent viscosity by 2 
orders of magnitude is able to significantly decrease the central rotation rate on the red-giant branch. 

Conclusions. Our results indicate that it is possible that the prescription for the horizontal turbulent viscosity largely underestimates 
its actual value or else a mechanism not included in current stellar models of low mass stars is needed to slow down the rotation in 
the radiative core of red-giant stars. 
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1. Introduction stars with significant convective envelopes (M+ < 1.4M ), the 
$H ' problem is complicated further by the magnetic braking by stel- 
es Stars rotate, and this rotation has important consequences on lar windfj rfKawalerlllQsk Another problem concerns the ini- 
their evolution. On the one hand, centrifugal acceleration re- tial rotation statej which can haye an influence on the rotation 
duces local gravity, mimicking a lower mass. On the otherhand^ profile during the pre . main seqU ence (PMS) and early main se- 
rotation induces meridional circulation ( Eddington|1925|; Mestd quence (MS) The MtM mmion state depends on the presence 
\1953 and shear and baroclimc instabilities ( |Mathis et al|2004|) , and nfetime of a circumste iiar disk during the first stages of the 
which contribute to the mixing of chemical elements. pMS; magn et i c fields seem to effective l y lock the stellar surface 
The problem of transport of angular momentum inside stars to the d j s k (IShu et alJll994h . iMaederf s l2009l textbook contains 
has not yet been fully understood. Several mechanisms seem to an exte nsive description of these processes, 
be active, the most commonly invoked being diffusion by tur- 
bule nt viscosity, transport by meridional circulation ( Maeder & Several stellar evolutionary codes already include transport 
Zahn ll998l:IZahnl|l992h. torques due to magnetic fields ("Maede r of angular momentum and associated chemical element mixing 
& Mevnet l2004l : lMathis & Zahnll2005l:IStrugarek et al.ll201lh . dChabover et al.lll995tlEggenberger et ~aT] l2008l: IF alacios et al l 
and transport by gravity waves dTalon & Charbonnelll2005h . For I2003I : iTalon et alj |1997l). Thev have been developed first with 

a focus on studying consequenc es on the evolution of massive 

Send offprint requests to: J.P. Marques stars or evolved stars (see, e. g., Maeder & Meynet 2000; Sills & 
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Pinsonneault 2000). The Sun and solar-like stars have been stud- 
ied as above (e. g JChabover et al Jl 1 995k iTalonll 1 997h IZahn et ail 
fl997h . and on e important application was the lithium surface de- 
pletio n (e. g.. ICharbonnel & Talonlll999l : lTalon & Charbonnell 
[l998h . 

One way to study the internal transport and evolution of an- 
gular momentum is to obtain seismic information on the internal 
rotation profile of stars in the low and intermediate mass range 
at different stages of evolu tion. The SoHO satellite for the Sun 
dTurck-Chieze et alj|2010h and the ultra high precisi on photom- 
etry (UHP) asteroseismic space missions CoRoT (Baglin 2003) 
and Kepler dBorucki et alj|2010l) offer such an opportunity. Seis- 
mi c information has being obtained for PM S and MS stars (e. 
g.. lBenomar et alJl20 lOt Escobar et al.ll2012h . and now we also 
have seismic observat ional constraints on the internal rotation of 
red-giant stars (e. g ., iBeck et all 120121 : iDeheuvels et al J 120121 : 
lMosseretai]l2012ah . 

Our goal is to use seismic diagnostics to test the descrip- 
tion of transport of angular momentum processes in 1-D stel- 
lar models of low-mass stars, from the PMS to the red-giant 
branch (RGB). Rotational splittings (differences between non- 
axisymmetric oscillation modes) are such diagnostics. When the 
rotation is slow, the relation between rotational splittings and ro- 
tation rate is linear and therefore relatively easy to interpret. 

The information that rotational splittings provide on the in- 
ternal rotation, however, depends on the physical nature of the 
stochastically excited modes, which in turn depend on the struc- 
ture of the star, hence its age. In addition, some rotation gradi- 
ent is expected to develop during the evolution; the core rota- 
tion rate is expected to become higher than the surface, as the 
core contracts and the envelope expands. If the rotation rate 
of some layers becomes too high, the associated distortion of 
the star and/or the Coriolis acceleration invalidates the linear ap- 
proximation that provides the usual (linear) rotational splittings. 
In such cases, non perturbative methods must be used (such as 
those developed i n iLignieres et al J 120061 : lOuazzani et alj|2012t 
iReese et al.l2006h . It is therefore necessary to determine if linear 
splitting are valid for an evolutionary stage in order to interpret 
correctly the observations. 

We have thus started a series of papers devoted to establish 
seismically validated processes of transport of angular momen- 
tum that play an essential role in shaping the rotation profile 
of stars. The present paper is the first of this series. The sec- 
ond and the third will focus specifically on the case of red-giant 
stars. They will investigate the seismic diagnostics of rotation for 
slowly rotating red-giant stars for which linear rotational split- 
tings are valid, and then rapidly rotating red-giant stars with non- 
perturbative methods. 

In the present paper, we compute the rotation profiles and 
their evolution with time from the PMS to the RGB with an evo- 
lutionary co de where rot ationally induc ed mixing using the pre- 
scription of lZahnl (1 19921) as refined by iMaeder & Zahnl dl998) 
has been implemented. We then follow the evolution of the lin- 
ear rotational splittings calculated with the rotation profiles ob- 
tained from our evolutionary models and discuss their validity. 

The paper is organized as follows: in Section [2] we describe 
the physical inputs that are implemented in our evolutionary 
code, with special emphasis on the transport of angular momen- 
tum and rotation induced mixi ng. We developed a version CE- 
SAM 2K of the code CESAM (lMorellll997t iMorel & Lebretonl 
2008) by implementing rotation-induced transport in radiative 
zones with careful attention to conservation and transport of an- 
gular momentum. The numerical scheme for the rotationally 
induced transport had to be modified compared to the general 



scheme used in CESAM2K, and some differences in the phys- 
ical inputs also exist. For these reasons, the modified version 
will be referred to as the CESTAM code hereafter. This will 
avoid any possible confusion with results from other versions of 
CESAM2K used in the international community. The acronym 
CESAM stands for Code d' Evolution Stellaire Adaptatif et Mod- 
ulaire, and the extra T in CESTAM stands for transport. A val- 
idated standard version of CESTAM will be made freely avail- 
able. 

The numerical implementation is described in Secfion[3] We 
then show comparisons with other codes for validation (Section 
|4). During the implementation of the code, we encountered some 
difficulties already present but not clearly explained in previous 
work. We chose to describe them carefully. In Sect. [5] we fol- 
low the evolution of theoretical linear rotational splittings for a 
1.3 M Q sequence of stellar models evolved from the PMS to the 
RGB. The ( = 1,2 modes are computed with the ADIPLS code 
dChristensen-Dalsgaardl l2008h in a frequency range where the 
modes are expected to be stochastically excited. We discuss the 
behavior of the splittings with frequency that can be observed 
during the various phases of evolution. For red-giant models, 
we find rotational splittings that are much higher than observed. 
In Section[6]we consider possible reasons for this disagreement. 
Finally, conclusions and some perspectives are given in Sect. [7] 



2. Physical input 



The original versio n of the CE SAM2K code dMorell 1997b Morel 
& Lebreton l2008l) consists of a set of routines that calculates 1-D 
quasi-hydrostatic stellar evolution including microscopic diffu- 
sion of chemical species. The solution of the quasi-static equi- 
librium is performed by a collocation method based on piece- 
wise polynomial approximations projected on a B-spline basis. 
For models without diffusion, the evolution of the chemical com- 
position is solved by stiffly stable schemes of orders up to four; 
in the convection zones mixing and evolution of chemicals are 
simultaneous. The solution of the diffusion equation employs 
the Galerkin finite elements scheme, and the mixing of chemi- 
cals in convective zones is then performed by a strong turbulent 
diffusion. 

The code CESAM2K allows the choice of several options 
for the physics. The microscopic input physics is updated regu- 
larly. The opacity tables are presently the OPAL95 data (Iglesias 
& Rogers Tl996l) c omplemented at low temperatures by the Wi- 
chita opacity data dFerguson et alJl2005h . Several sets of opacity 
tables are provided that corres pond to various mix t ures o f chem- 
ical elements, for instance the Grevesse & Noels ( 1993) or As- 
plund et al. (120051 12009) solar mixtures, as well as a-element 
enhanced mixtures. We included P othekhin's updated conduc- 
tive opacities (see e. g. ICassisi et al.ll2007l) . Several options 
are possible for the equation of state, the most comm only used 
being the OP AL2005 EoS dRogers & Navfonovl 120021: Rogers 
et al. 119961) . Several networks of nuclear reactions (and corre- 
sponding updated nuclear reaction rates) were implemented, al- 
lowing the evolution of stars to be calculated from the PMS up 
to helium burning. The microscopic diffusion of chemical ele- 
ments includes gravitational settling, thermal, and concentration 
diffusion terms, but no radiative accelerations. T wo formalisms 
are available for microsc opic diffusion transport (Burgers 1969; 
iMichaud & ProffitS 1 1 993h . Also, two options a re available for 
treatin g conv ection, the classical MET theory of Bo hm-Vitensd 
d!958l) or the lCanuto et alj d 19961) so-called full spectrum of tur- 
bulence. The atmospheric boundary condition is derived either 



Article number, page 2 of 1 151 



J. P. Marques et al.: Seismic diagnostics for transport of angular momentum in stars 



from gray mod el atmospheres or from Kurucz ATLAS9 models 
(Kurucz 2005). Overshooting is an option. 

Mass loss can be considered with different prescriptions. In 
the mod els presented he reafter we use the empirical mass loss 
rates of Reimersl (| 19751) scaled with metallicity according to 
ISchaUeretalJ(ll992h . 

Recently, the CESAM2K code has been involved in the 
ESTA activities undertaken to prepare the interpretation of the 
seismic observations of CoRoT. Stellar models were calculated 
for a range of mass, chemical composition, and the evolution- 
ary stages corresponding to CoRoT main targets and have been 
compared with the results of several o ther evolutionary cod es 
showing a very good g eneral agreement dLebreton et al J2008allr3 : 
iMontalban etafll2008h . 



CE SAM2K is freely ava ilable for download, with the details 
in lMorel & Lebretonl(l2008l) . 



2. 1 . Transport of angular momentum 

In convective zones, although there is differential rotation in lat- 
itude, the mean rotation rate at a given radius weakly depends 
on the radius. Therefore, it is often assumed in 1-D stellar evo- 
lution codes that convective zones rotate as solid bodies (e. g . 
iMevnet & MaedeifcOOOtlPalacios et ai1l2003l:lTalon et al.lll997h . 
In extended convective zones, however, comparisons with 3-D 
numerical simulations suggest inste ad a prescription of uniform 
specific angular mo mentum (e. g. Deniss enkov & Toutl feOOO: 
IPalacios et ai1l2006h . Both options can be used in CESTAM. 

In a r adiative zone, we used the formalism of Zahn ( 1992), 
refined in Ma eder & Zahnl dl998h . to model the transport of an- 
gular momentum and chemical species by meridional circulation 
and shear-induced turbulence. In what follows, we sketch the 
model to make it clear exactly which equations we adopted, as 
several slightly different versions have been used (e. g. Denis- 
senkoy & Tout 120001: iMaeder & Zahnlll998HMathis & Zahnl 
I2004t IPalacios et al.ll2003HTalon et al.ll 1 997t IZahnll 1 992b . 



Turbulence is expected to be highly anisotropic owing to the 
stable stratification in radiative zones. Turbulence would then 
be much stronger in the horizontal than in the vertical direction. 
Thus, the hypothesis of "shellular rotation" can be used: as dif- 
ferential rotation is presumably weak along isobars, it is treated 
as a perturbation. All variables / can b e split into a mean value 
over an isobar and a perturbation (as in Zahn 1992): 



/(p,0) = /(p)+/a(p)Pz(cos0), 



(1) 



where P2(cos §) is the second-order Legendre polynomial and p 
is the pressur e. Higher order effect s, included in the formalism 
developed by Mathis & Zahn ( 2004), are not considered here for 
the moment. 

The velocity of meridional circulation can be written in a 
spherical coordinate system as 



-> _> diVcos #)_> 

U = U 2 (r)P 2 (cQsffrfr + V 2 (r) 2 \ e, % 

d# 



(2) 



where r is the mean radius of the isobar and § the colatitude. 
The vertical compone nt U 2 is given in Appendix [A] following 
iMaeder & Zahnl d 1 998h . The equation of continuity in the anelas- 
tic approximation gives then the horizontal component V 2 '. 



where p is the density. 



(3) 



The transport of angular momentum obeys an advection- 
diffusion equation: 



P 



dt 



— - pr 4 OC/ 2 ) + -- 
5r z or v 1 r i Q- r 



4 dCl\ 



dry 



(4) 



where vy is the vertical component of the turbulent viscosity and 
d/dt represents the Lagrangian time derivative. 

The relative horizontal variation of the density, = p/p, and 
the mean molecular weight, A = pip, obey 



= 
dA 



2r2 Q<9 Q 
3g dr 



U 2 



6D h 



Hp 



A, 



(5) 
(6) 



where is the horizontal component of turbulent diffusion dis- 
cussed in Appendix [B] The evolution of A, Eq. ©, depends 
on the competition between the advection of a mean molecular 
weight gradient and its destruction by the horizontal turbulent 
diffusion. 



2.2. Evolution of the chemical composition 

The vertical advection of chemicals due to the large-scale merid- 
ional circulation coupled with a strong horizontal turbulent dif- 
fusion results in a vertical diffusion process (see, e. g. Chaboyer 
& Zahn 1 19921) . The equation of the chemical composition evo- 
lution can then be written as 



d£i = 8_ 
dt dm 



(4nr 2 pf (Dy + D eff ) 



dXi 
~dt 



dXi 

nucl dt 



(7) 



where X-, is the abundance by mass of the z'-th nuclear species 
and Dy = vy and D e ff the vertical diffusivity and the diffusion 
coefficient associated with the meridional circulation; Dy and 
D e ff are discussed in AppendixlBl 

A necessary con dition for shear instability is the Richardson 
criterion as given bv lTalon & Zahnl d 19971) . Another condition is 
that the turbulent viscosity vy must be greater than the molecular 
viscosity, v, as expressed by the Reynolds criterion 



vv > vRe c , , 

whe re Re,- = 
et al.T2000). 

Vy = V. 



(8) 



: 10 is the critical Reynolds number (Schatzman 
When condition ® is not satisfied we use Dy = 



2.3. Initial conditions 

Stars are fully convective when they start their PMS evolution 
on the Hayashi track. Assuming that convective zones rotate 
like solid bodies, the star should have uniform angular velocity 
at the beginning of its evolution. 

Several facts complicate this simple picture. First, Palla 
& Stahler (119911) showed that stars that are more massive than 
about 2M Q are no longer fully convective when they appear on 
the PMS (the exact mass depends mainly on the protostellar ac- 
cretion rate and deuterium abundance). Second, if during the 
PMS the only process slowing down stellar rotation were the 
magnetic braking by stellar winds mentioned above, stars would 
rotate much more rapidly than observed on the ZAMS. The PMS 
is too short for this process to slow down the star significantly. 
An additional process is needed du ring the PMS, most likely disk 
locking (e. g jBouvier et al.l[l997h . 
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Young stars are most often surrounded by a circumstellar 
disk left over after the main accretion phase is over. The mag- 
netic coupl ing between the star and the disk slows the star down 
(see, e. g.. lShu et al.lll994l) . and this effect is often modeled by 
assuming that the stellar surfa ce corotates with th e disk at a con- 
stant angular velocity (see iBouvier et al. I fl997h as long as the 
disk exists. Once the disk has disappeared, the star surface ro- 
tation evolves freely. This scenario is implemented in CESTAM 
with the disk lifetime and period fdisk as free parameters of 
the model. 

2.4. Magnetic braking 

Stars that are less massive than about 1.4M have significant 
outer convective zones. A solar-type dynamo operates there and 
generates a magnetic field. The coupling between the magnetic 
field and the plasm a in the stellar wind strongly brakes the ro- 
tation of the star. iKawalerl (Q988, hereafter K88) proposed the 
following law for the loss of angular momentum /: 

l _i 

J = -K w QOL(|-) 2 (|-) 2 forQ>Q sat , (9) 

where £2 sat is a saturation angular velocity, above which mag- 
netic field generation seems to saturate. The v alue of fi sat is 
often set at 8 - 14Q (as in IBouvier et al. I fl9971) . There are in- 
dications, however, that f2 sat varies with stellar m ass, (see, e. g. 
lAndronov et alJl2003t iKri shnamurt hi et aill 19971) . 

The parameter in Eq. (0 is usually calibrated by re- 
quiring that calibrated solar models have Q = £2 = 2.86 x 
10" 6 rads~'. The precise value of needed to spin down the 
Sun to its current period depends on the prescription for the 
transport of angular momentum adopted (see section l4~2l below). 
The parameter A"w should also depend on stellar mass (see, e. g. 
lKrishnamurihietal.il 19971). 

Recently. iReiners & Mohantvl d2012l hereafter RM12) have 
criticized the approach used in K88. Specifically, K88 supposed 
that the surface magnetic flux goes as some power of the angular 
velocity, whereas RM12 suggest instead that it is the magnetic 
field strength that obeys such a law. As a consequence, they 
obtain the following law: 



J = 3 forO<Q sat 

j = -KmsQCtt^j 3 l^j 3 for£2>Q sat . (10) 

The value of Krm does not depend on the stellar mass, and Q sat 
3Q . Both prescriptions, Eqs. (0 and ([Tol l, are implemented in 
CESTAM and tested below. 

3. Numerical procedure 

3.1. The overall problem 

Stellar evolution depends on two interconnected problems: the 
problem of stellar structure (solving the stellar structure equa- 
tions), and the problem of the evolution of the chemical compo- 
sition. It has proved very difficult to solve the two problems 



simultaneously; different stellar evolution codes employ sev- 
eral techniques to overcom e this difficulty. Some codes (e. g., 
iDeglTnnocenti et alj|2008l) compute the solutions of the struc- 
ture equations and chemical composition in an independent way, 
where they use the structure calculated at the previous time step 
to evolve the chemical composition to the current time step, and 
then use this composition to calculate the structure at the cur- 
rent time step (or the other way around: first the structure, then 
the chemical composition). A second kind of code computes 
the chemical co mposition between eac h iteration of the structure 
probl em, as in Scuflaire et al. | d2008l) . And finally, Eggleton's 
code (Eggleton| [l971l) solves the two problems simultaneously. 
IStanc liffe (2006) has shown that differences between the three 
kinds of codes are only significant on the AGB. 

The CESAM2K code dMorelll997l:lMorel & Lebretonl2008l) 
belongs to the second category above. It begins a time step 
by updating the chemical composition. With the new chemi- 
cal composition, a new structure is obtained by performing one 
iteration of the algorithm used to solve the structure equations. 
This structure is used to update the chemical composition again, 
before a new iteration on the structure is performed. The pro- 
cedure is repeated until convergence of the structure algorithm. 
CESTAM keeps this iterative scheme. 

Rotation with angular momentum transport introduces a new 
problem interconnected with the previous two. The evolution 
of the chemical composition depends on the turbulent mixing 
induced by differential rotation, whereas the structure equations 
are changed by the inclusion of the centrifugal acceleration. In 
CESTAM, we inserted the resolution of the angular momentum 
transport equations at the beginning of the cycle, because the 
rotation profile is required to calculate the turbulent diffusion 
coefficients that are needed to update the chemical composition. 
Our cycle, then, is as follows: we update the rotation profile, then 
the chemical composition, and finally we iterate on the structure. 
The cycle is repeated until convergence. 

In the absence of external torques, total angular momentum 
is conserved. The rotation profile obtained at this stage does not 
enforce angular momentum conservation, however, because the 
profile O(m) was computed with the stellar structure obtained in 
the iteration before convergence. The total angular momentum 
is / oc f Q,(m)r 2 (m)dm and the function r(m) changes between 
iterations. We need to compute the rotation profile one more 
time after convergence to make sure that it is consistent with 
the stellar structure, so that angular momentum is numerically 
conserved. 

3.2. The rotation profile 

Equation (01 (with Eq. dA.ll ). Eq. © and Eq. ©X is a fourth- 
order differential equation in r. To solve it, we split it into four 
first-order differential equations , and solve the system using the 
well k nown relaxation method dHenvev et alJI 19641 : iPress et ail 
2007). Equation (0 is solved simultaneously with Eq. (0), so 
that we have a total of five finite difference equations. We chose 
this method for its simplicity in dealing with the complex rela- 
tions between variables expressed in Eq. (0), taking Eqs. dA.ll ). 
(0), and © into account. Methods that require that the solutions 
are approximated by a linear combination of known functions 
(collocation methods, spectral methods) are therefore difficult 
to implement. In practice, the relaxation method we employed 
proved to be efficient, robust, and fairly stable. 

The scheme can be fully implicit or semi-implicit in time 
(to ensure higher order accuracy in time). The time step can be 
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subdivided if needed, a useful feature for future developments 
involving faster processes. 

Four boundary conditions are needed. At the top of a radia- 
tive zone, we impose conservation of angular momentum and no 
differential rotation: 



Table 1. Parameters of calibrated solar models. 



^=0 
dr 



II >) 



Am 



4n , 
—pr 4 D.U 2 



(ID 



where m, is the mass inside the top of the radiative zone. If the 
convective zone above it is at the surface, j is the torque applied 
at the surface of the star, Eq. (0, otherwise j = 0. At the bottom 
of a radiative zone, similarly, 

dr 

- (r 2 Q)dm = — pr A Q.Uo , (12) 
dfJo V 1 $ m=m b 

where mj is the mass inside the bottom of the radiative zone. If 
the center of the star is radiative (m/, = 0), Eq. (TT2l) is replaced 
by Ui = 0, the requirement that there is no mass flow out of the 
center. 

Intermediate convective regions (regions that are neither at 
the center nor at the surface) are treated as if they were special 
points. The equations between the beginning and the end of an 
intermediate convective zone (between m\ and mi) are 

d /"* mf 47T 

-J (r 2 ti)dm = y^ifn^y-Pf^ntUv) ( 13 ) 



<9Q 



= 



(14) 



dQ. 
~dr~ 


f 


= 


(15) 


A 




= 


(16) 


A f = 


= 0. 


(17) 



Equation ( 113b guarantees conservation of total angular momen- 
tum, Eqs. (1141115b impose no shear at the borders of the con- 
vective zone, and Eqs. (1161117b result from the absence of jj- 
gradients in a convective zone. 

The four equations resulting from Eq. (|4]l are written in 
finite-difference form between pairs of points, say, aim - m.k-\ 
and m = m\. Equation (0 requires a special treatment, because 
it is not a differential equation in r (or m). We write Eq. © as 
an algebraic equation at point m = m&_i, because otherwise we 
would be solving for (At + Aj_i) /2 (the average between m^\ 
and mil instead of the values at mk-\ and m^). Indeed, A only 
appears as (A* + A^-i) /2 in the other equations, since it has no 
space derivatives. A fifth boundary condition is needed: A = 
at m — m, (the top of the radiative zone) as in Eqs. ( 1161117b . 

4. Comparison with results from other evolutionary 
codes 

The complicated nature of the equations describing the evolution 
of the rotation profile (and particularly the evolution of Ui) made 
it important to validate our approach. To do that we compared 
our results with those obtained with other implementations of 
rotational mixing in s tellar evolution codes, namely STAREVOL 
(Palacios et al. 2003) for the solar case (see also Turck-Chieze 
et al. 2010) and the Geneva stellar evolution code (iTalon et al.1 
1 19971) for the case of 3M and 5M stars. 



Case 


Q'CGM 


Yi 


Zi 


K w (cgs) 


GN93, no b. 


0.7031 


0.2792 


0.02101 





GN93, b. 


0.6842 


0.2760 


0.02019 


6.2 x 10 47 


AGS05, no b. 


0.6702 


0.2607 


0.01371 





AGS05, b. 


0.6502 


0.2566 


0.01307 


5.6 x 10 47 



Notes, q-cgm is the mixing length, Y[, Z\ the initial helium and metallic- 
ity abundances, and A'w the braking parameter entering Eq. l[9}. No b. 
indicates no braking, b. braking. 

Table 2. Characteristics of calibrated solar models. 



Case 


Y s 


z s 


rcz/Re 


GN93, no braking 


0.2511 


0.01950 


0.7123 


GN93, braking 


0.2613 


0.01950 


0.7166 


AGS05, no braking 


0.2319 


0.01260 


0.7314 


AGS05, braking 


0.2425 


0.01260 


0.7364 



Notes. Y s and Z s are the surface helium abundance and metallicity, re- 
spectively, and r C zlRe ' s me normalized radius at the base of the con- 
vective envelope. 



4.1. Standard physics 

In these comparison s, we used the OPAL equation of stat e 
dRogers et al.lll996l) and opacities (Iglesias & Rogers 1996), 
compl emented at T < 10 4 K by the lAlexander & Ferguson! 
(1994) opacit ies. W e used the NACRE nuclear reaction rates of 
Angu lo et al.l ([T999) except for t he 14 N+p reaction, for which we 
used the reaction ra tes given in llmbriani et all d2004h . The so- 
lar compositions of iGrevesse & Noelsl (I1993L hereafter GN93) 
and Asolund et all (I2005L hereafterAGS05) were used. The 
Schwarzschild criterion was used to determine convective insta- 
bility. Convective core overshoot fully mixes the chemical com- 
position to a distance d m - a m x min(Hp, r co ) from the border 
of the convective core (where r co is the radius of the core deter- 
mined by the Schwarzschild criterion). The temperature gradient 
in the overshoot zone is V = V a d. The centrifugal acceleration 
is taken into account by adding the average centrifugal acceler- 
ation 2f2 2 r/3 to gravity in the hydrostatic equilibrium equation. 
The atmosphere is computed in the gray approximation and in- 
tegrated up to an optical depth of r = 10~ 4 . 



4.2. The Sun 

We computed calibrated rotating solar models using the compo- 
sitions of GN93 and AGS05, with and without magnetic brak- 
ing according to the magnetic braking law, Eq. {9). Models 
were calibrated to within 10~ 5 in luminosity, radius, and surface 
metallicity (Z = 0.0195 for GN93, Z = 0.0126 for AGS05). Ini- 
tial rotational velocities (and parameter Kyj for the models with 
magnetic braking) were chosen so that models have an equatorial 
velocity u eq = 2.02 kms" 1 at the solar age (4.6 Gyrs). 

The temperature gradient in convective zones was computed 
using the convection model of Canu to et al.1 (Q996) with I — 
&cgmHp. All models include microscopi c diffusion and set- 
tling u sing the approximations proposed by Michaud & Proffitt 
(Il993l) . The initial rotational velocities for the cases with brak- 
ing were chosen in order to have u eq = 20 km s _1 at the ZAMS. 

Table[T]shows the values of parameters resulting from the so- 
lar calibration in ccgm, the initial helium, and metallicity abun- 
dances (Y[ and Zi). The greatest differences between the cases 
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Fi e. L Differences in soun d speed profiles between the seismic models 
of iTurck-Chieze et al.1 feOOlj) and CESTAM. 



with and without magnetic braking concern the initial abun- 
dances. This is because mod els with braking rotate m uch faster 
at the center (as shown in ITurck-Chieze et al.ll2Q10h and thus 
have a higher Q-gradient, leading to a much higher turbulent 
diffusion. Turbulent diffusion tends to homogenize the radiative 
zone to a large extent, partially erasing the composition gradient 
created by gravitational settling. 

Table [2] shows some characteristics of the calibrated solar 
models. The higher turbulent diffusion partially stops the settling 
of helium, causing a higher helium abundance in the convective 
zone. However, the strong O-gradient predicted by models with 
braking in the radiative zone is in direct contradiction with he- 
lioseismic results, which indicate a flat rotation profile to within 
r - 0.25Rn of the center (see, e. g. Chaplin et al. 1999; Couvidat 
et al. 20031 : lEff-Darwich et al.l2008l) . A newphvsical mechanism 
is needed to explain the discrepancy, such a s transport of angular 
mom entum by internal gravity waves (see Talon & Charbonnel 
l2005h a nd/or magnetic stresses dMaeder & Mevnet)2004t Mathis 
&Zahn l2005T ~ 

Figure Q] shows differe nces between the sound spe ed c s pro- 
files of the seismic Sun of ITurck-Chieze et al ] (120011) using re- 
sults from SoHO (GOLF-MDI) and our calibrated solar models 
(see also lGoupil et al. 1 1201 1| for intermediate o pacities AGS 09). 
The seismic Sun of Turck-Chieze et al. (2001) agrees with Basu 
et al. (2009) using results from BiSON and MDI. A good agree- 
ment at the level of 0.3 % is seen when using the old mixture of 
GN93 with only small discrepancies below the convection zone 
and in the central region. On the other hand, severe discrepan- 
cies occur when more recent mixtures are used, such as AGS05. 
These resu lts are quite similar to what is found i n the literature 
(see, e. g.. iBasu & Antiall2008HBasu et alJll997t Turck-Chieze 
et al. l201 ll for models with no rotation induced transport). 

We compared our results for calibrated solar models includ- 
in g rotation-induced transport of Type I with the results obtained 
in ITurck-Chieze et al.1 d2010l) with the STAREVOL code. We 
found good agreement between the profiles of the angular ve- 
locity, U2, and the two components of turbulent viscosity. Sev- 
eral improvements have been implemented i n the code CESTAM 
compared to the CESAM2K version used in lTurck-Chieze et al.1 




Fig. 2. Profiles of the vertical component of the meridional circu- 
lation E/a(r) as a function of normalized radius (r/i? s ) for 5 M Q stel- 
lar models computed with CESTAM (continuous line) and the Geneva 
code (dashed line) when the central hydrogen abundance is X c = 0.35 
(see text). The left panel shows the central regions, the right panel the 
surface. 



(120101) . but with no consequence for the solar case. We confirm 
that rotation induced transpor t of Type I does no t help remove 
the discrepancies, as found bv lTurck-Chieze et al.l (120101) . 



4.3. Higher mass, main sequence stellar models 

We compared results for models with higher mass stars. These 
stars have thin convective or fully radiative surface layers and 
are therefore not expected to undergo magnetic braking. Cal- 
culations here are carried out assuming global conservation of 
angular momentum, which is obtained at the precision level of 
10~ 6 . 

The first comparison concerns the evolution of a 5M H model 
comp uted with CESTAM and Geneva codes (as in iTalon et akl 
1 1997b . using X{ = 0.73 and Z; = 0.01 (and no magnetic braking). 
These models were computed without microscopic diffusion and 
settling. Figure [2] shows the profile of U<i(f) at the middle of the 
MS (when X c = 0.35). Curves for the same quantities superim- 
pose, showing excellent agreement between the results obtained 
with the two codes. There are small differences that we attribute 
to the different microphysics used in the codes. The rotation pro- 
files, not shown, also agree during the course of the evolution. 

We also compared the effects of rotational mixing on the evo- 
lutionary tracks on the Herzsprung-Ru ssel diagram (HRD) be- 
tween CESTAM models and those of Eggenber ger et al.| (120101) 
computed with the Geneva code dEggenberger et al. 2008). Fig- 
ure [3] shows evolutionary tracks on the HRD for a 3M model 
calculated without extra mixing (microscopic diffusion, convec- 
tive overshoot, rotation), with rotational mixing only and with 
overshooting only. 

We considered two cases for the models with overshooting, 
a ov = 0.1 and a m = 0.2. The model with rotational mixing ha s 
u eq = 150 km s" 1 attheZAMS. As in lEggenberger et al.l (12010). 
the evolution of the rotating model closely resembles the evolu- 
tion with a m =0.1. The difference between the rotating models 
and the nonrotating models at the ZAMS is due to the centrifugal 
acceleration. The decrease in the effective gravity (gravity mi- 
nus centrifugal force) replicates a nonrotating star with a lower 
mass. 
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Fig. 3. Evolutionary tracks on the HRD for a 3M Q model calcu- 
lated with CESTAM code assuming either no extra mixing (dashed) 
or rotation with an initial equatorial velocity u cq = 150kms~' at the 
ZAMS (continuous line) or without rotation but with overshooting and 
a m = 0.1 (dot-dashed line) and a m = 0.2 (dotted line). 
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Fig. 4. Evolutionary track for a 1.3M Q sequence of models from 
PMS to RGB. Dots indicate models for which rotational splittings are 
computed. The numbering of the selected models stars with the first dot 
on the PMS and ends with number 25 as the highest dot on the red-giant 
branch 



The ev olutionary tracks shown in Fig. [3]are identical to those 
in Fig. 2 of Eggenberge r et al.l d2010h . 



4.4. Evolved low-mass models 



The rotation profiles that we obtain for our red-g i ant m odels are 
similar to those computed by Egg enberger et al.l d2010l) . Exam- 
ples for M+ = 1 3M Q are given in the section below. 



5. Evolution of rotational splittings from PMS to 
RGB 

Rotational splittings are useful as seismic diagnostics for mea- 
suring the stellar internal rotational profile. For slow rotation, 
a first order perturbation description provides rotational split- 
tings that are linearly dependent on the rotation profile. They are 
therefore convenient tools, easy to compute from stellar mod- 
els and theoretical oscillation codes for comparison with ob- 
servations. The first order approximation provides the follow- 
ing expression for the linear rotational splittings (Christensen- 
Dalsgaard & Berthomieu ll991 , and references therein): 



5v„ e = I K„, 
Jo 



e(x) -z — Ax 
2n 



(18) 



where x = r/R* is the normalized radius, O is the angular rota- 
tion (rad/s) and the rotational kernel K„e takes the form 



Knt = ^-\fr+ W + l)? h ~ 2{S, - 

where I„i is the mode inertia 

he = f [f r + £(£ + l)f h ] P x 2 Ax. 



p x 



(19) 



(20) 



The quantities entering the equations above are the fluid vertical 
and horizontal displacement eigenfunctions, £ r and ^ respec- 
tively, and the density p. For solid-body rotation, the rotational 
splittings become: 



6v n e = p n c 



with 



2^ 



(21) 



1 f 1 

pV = 1 -C nt = — \f r +{({+ l)f h -2^ h -f 2 ] p x 2 Ax, (22) 

int JO 

where C„t are the Ledoux coefficients. 

For asymptotic pure p-modes, C n ( ~ and fi n ( ~ 1 . For pure 
g-modes,C„f ~ l/i(i+l)anA/5 nC ~ 1 - 1 /€(£+ 1), i. e. J3 n£ ~ 1/2 
for £ = 1 modes. 

To investigate the evolution of rotational splittings with stel- 
lar age, we computed a 1.3 M evolutionary sequence including 
rotationally induced transport as described in Section [2] above. 
We used the same input physics as in the case AGS05 with brak- 
ing described in Section \4~2\ but without microscopic diffusion 
and settling, as we found that microscopic diffusion and set- 
tling do not affect the rotation profiles of our models. We chose 
M* = 1.3M Q because it is a typical seismic mass obtained for 
the Kepler and CoRoT red-giant stars. 

Figure |4] shows an evolutionary track for a sequence of 
1 .3M stellar models in a HR diagram computed with CESTAM. 
In order to interpret the rotational splittings, we plot the evolu- 
tion of the central and surface rotation rates with age in Fig. [5] 

We then computed oscillation frequencies and rotational 
splittings for several models spanning the track (shown in Fig. |4j 
using the freely available AD IPLS adiabatic oscillation code 
(Chri stensen-Dalsgaardl [2008h . We computed frequencies of 
£ — 0, 1 axisymmetric modes and rotational splittings for £ = 1 
modes. 

The 1 3M Q models have a convective envelope from PMS to 
the RGB, hence we expect stochastically excited modes all along 
the sequence. Thus, we chose the frequency range that spans an 
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Fig. 5. Evolution of the central (blue line) and surface (green line) 
rotation rates for a 1.3M Q sequence of models from the PMS to the 
RGB. The insert shows the central and surface rotation rates along the 
red-giant branch. 



interval of a few radial orders below and above n max , the radial 
order corresponding to the frequency at maximum power. It is 
estimated as n max = v max /Av, where the frequency at maximum 
power spectrum v max and the mean large separation Av are given 
by the usual scalings relations 



Av 



M ' K R 

v l/2 



eft' 



v-1/2 



eif,0 



(23) 



(24) 



with r eff , G = 5777/T, v max = 3050/rfIz and Av G = 134.7//HZ 
for the Sun. Indeed, it has been conjectured, then shown obser- 
vationally, that these relations predict well the location of the 
excited frequency range of stochastically excited modes (e g. 
Brown et al. 1 11991b iKallinger et ail 1201 Ob iKieldsen & Beddind 



19951) 



5.1. Validity of linear approximation for the rotation splittings 

The validity of the linear approximation is estimated by compar- 
ing the rotation rate to the oscillation frequency. In a perturbative 
approach, the parameter 



2D. 



(25) 



is thus assumed smaller than unity. We evaluate this parameter at 
the center of each selected equilibrium model of the evolution- 
ary sequence, where Q is largest, and for a range of frequency 
spanning the radial order interval (n max - 4, n max + 4). 

Sharp rotation gradients develop in the central layers and can 
also cause some departure from the linear approximat ion. In the 
oscill ation equations (see, e. g. Eqs. (34.7-34.12) in luiino et al.l 
1989), the rotation gradient term appears with a factor 



1 



dQ 2 



(oj + nfdlnr' 

Thus, we consider that the linear approximation is valid if 
1 dQ 2 1 dQ 1 dQ 



(26) 



(cj ± CI) d mr w + Qdlnr w ± Q d In r 



« 1, 



(27) 
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Fig. 6. The evolution of the parameters f, Eq. J25t (top) and A, 
Eq. J28t (bottom), with the effective temperature of the model along its 
evolution from the PMS to the red-giant branch. To be conservative, the 
layer where CI is maximum is considered for f and the layer where the 
gradient is maximum is chosen for computing A. For (, we used m = 
-1, since £ is higher for prograde modes, while for A we considered 
m = 1 as no difference is found for m = - 1 due to the fact that CI « o> 
where A is maximum. 



where the first inequality follows from assuming £ < 1 . 
We then evaluate the quantity 



_ / dQ. x 1 
= Ulnr^ZT+Q 



(28) 



at the radius where the rotation gradient is greatest and for 
<x>l2n = Vmax. Fig. [6] shows both quantities, f and A, as func- 
tions of the effective temperature of the models. Both parame- 
ters remain much smaller than unity until the model reaches the 
base of the red-giant branch. Therefore, the use of the linear ap- 
proximation for the rotational splittings can be safely used for 
stochastically excited solar-type modes of low mass stars except 
for the fast rotating cores of red giants. 

5.2. Evolution of rotational splittings of stochastically excited 
modes along an evolutionary track 

In the following, we discuss the information that can be retrieved 
from the average value throughout the star (Q/2tt) = Sv„f//3„e 
weighted by the mode. This definition is given at fixed t, and 
depends on the radial order n. 

- The PMS regime: On the PMS, the central and surface rota- 
tion rates are equal as long as the model is fully convective. 
The rotation rate remains constant in time for the duration 
of disk locking. When the radiative core appears, the central 
rotation starts to increase due to the contraction of the cen- 
tral layers. When disk locking stops, the surface rotation rate 
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Fig. 7. Top panel: Rotation profiles as a function of the normal- 
ized mass m/Mi, for l.3M e PMS models # 1 (black), # 2 (cyan), # 3 
(red), # 4 (green), # 5 (blue) and #6 (magenta) shown in Fig. [4] Sec- 
ond panel: Rotational splittings for t = 1 modes as a function of the 
normalized frequency v/Av for models # 1,# 3,# 5 and #6 (same col- 
ors as above). Third panel: B„c for the same models. Bottom panel: 
{Q./2n) = Sv„{/B„ e for the same models. The large separation Av goes 
from 87.9 pHz (model # 1) to 106.2 pMz (model # 6). 



Fig. 8. Top panel: Rotation profiles as a function of the normalized 
mass m/M* for 1.3M G MS models # 6 (magenta), # 7 (blue), # 9 (red), 
#10 (green) and #11 (black) shown in Fig. [4] Second panel: Rotational 
splittings for E = 1 modes as a function of the normalized frequency 
v/Av for models # 7,# 9,# 10 and #11 (same colors as above). Third 
panel: B„[ for the same models. Bottom panel: {Q.) = Sv„c/B„( for the 
same models. The large separation Av goes from 97.6yuHz (model # 7) 
to 63.2/iHz (model* 11). 



starts increasing as the star contracts. Fig. [7] shows the ro- 
tation profiles of the selected PMS models shown in Fig. [4] 
The surface rotation rate evolves from 2.84 fiHz (model # 
1) to 3.9 //Hz (model # 5) while the central rotation rate in- 
creases from 1 1.3 //Hz to 21.2 //Hz. The rotational splittings 
for ( — 1 modes are computed for models # 1 to # 6 accord- 
ing to Eq. (118b and are also shown in Fig. [7] The splittings 
are nearly independent on radial order and increase steadily 
with the surface rotation rate as the model evolves along the 
PMS towards the ZAMS. The value of the Ledoux constant 
almost vanishes, and j3„c ~ 1, because the modes are essen- 
tially pressure modes in the frequency range where they are 
expected to be stochastically excited. These modes have no 
amplitudes in the inner layers and cannot probe the rotation 
there. Thus the mean rotation rate (Q.)/2n corresponds es- 
sentially to the rotation averaged over the outer layers. As 
the rotation increases inwards in these layers, we find that 
{Q)l2n = 3.8 yuHz for Model # 1, close to - but slightly larger 



than - its surface rotation rate. Model # 6 is already on the 
MS. Its surface rotation is faster than that of younger mod- 
els but its central rotation is slower. This is not seen in the 
corresponding splittings and the averaged rotation {Q.)/2n 
remains close to the surface value. Results are similar for 
{ = 2 modes (not shown). 

- The MS regime: On the MS, the surface rotation rate de- 
creases with age due to both braking at the surface and an 
increase of the stellar radius. The central rotation rate also 
decreases with time due to transport of angular momentum 
from the core to the surface. The rotation in the central re- 
gions first decreases from Model #6 to #9 and increases by 
roughly 60 % from model #9 to #1 1 (Fig. [8). The excited 
modes are still in the frequency domain of high order p- 
modes and their rotational splittings reflect the surface be- 
havior only. The splittings decrease with the surface angu- 
lar velocity. At the end of the MS (model #11), the central 
rotation rate has already reached 15.4 fiHz, while the sur- 
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Fig. 9. Top panel: Rotation profiles (logarithmic scale) as a func- 
tion of the normalized mass m/M for 1.3M Q subgiant models # 12 (blue 
solid) to # 17 (black solid) as indicated in Fig. [4] Second panel: Ro- 
tational splittings for I = 1 modes as a function of the normalized fre- 
quency v/Av for model #12,#14,#15 and #16. Third panel: B„ ( for the 
same models. Bottom panel: {Q.)/2n = Sv„c/B„t for the same models. 
The large separation Av goes from 56.5 yuHz (model # 12) to 38.0//Hz 
(model #16). 



face rotates at a rate of 0.6 /iHz: the core is rotating roughly 
26 times faster than the surface. As can be seen in Fig. [8] 
the rotational splittings for this model yield a mean rotation 
(Q)/2tt = 3.8 /iHz, that is 1.22 /iHz higher than the surface 
rotation rate but much lower than the central rotation rate. 

The subgiant regime: When the model cools along the sub- 
giant branch, its rotation keeps evolving differently in the 
inner and outer parts, slowing down at the surface and accel- 
erating in the inner regions (Fig. [9). Substantial changes 
appear when the model leaves the MS and evolves as a 
subgiant. The Brunt-Vaisala frequency of red-giant stars 
is very high in the radiative interior, and the frequencies 
of g-modes enter the frequency domain of p-modes in the 
frequen cy range where stochastically excited modes are ex- 
pected dChristensen-Dalsgaardll201 lri . During the subgiant 
phase, avoided crossings between p- and g-modes appear. As 
observ ationally shown dDeheuvels et alj|2010t : iMosser et alj 
2012b), some p-modes become mixed modes. These modes 



Fig. 10. Rotation profiles as a function of the normalized mass m/M^ 
for 1.3M Q models # 18 to # 21 and # 23, # 25 (from blue to black) on 
the red giant branch as indicated in Fig. [4] 



have amplitudes both in the inner regions, where they behave 
as gravity modes, and i n the surface layers, where they be - 
have as acoustic modes dDziembowski & PamvatnvkhTl 991). 
They are very interesting as they allow probing regions deep 
inside the star. The number of such mixed modes is small 
at the beginning of the subgiant branch and increases as the 
star evolves toward the red-giant branch. During this phase, 
modes will eventually undergo avoided crossings, exchang- 
ing nature from p to g. 

Figure|9]shows that Model #12 is not evolved enough so that 
no mixed modes are present in the frequency domain shown. 
The splittings remain nearly frequency independent and pro- 
vide a mean rotation rate close to its surface value. The vari- 
ation in the splittings with frequency for the more evolved 
models # 13 to # 17 clearly shows the presence of mixed 
modes. With amplitudes in the dense interior, mixed modes 
have a much stron ger inertia than the neighboring p-modes 
(Dupr etet al.l20 09). For the same reason, the rotational split- 
tings exhibit the same behavior as the mode inertia although 
in a more pronounced way. They vary from one mode to the 
next. They are also larger than for pure p-modes when the ro- 
tation increases inward within the model. For subgiant stars, 
only one I = 1 mixed mode exists between two consecutive 
radial modes. This explains the saw-like aspect of the varia- 
tion in the splittings with frequency. This is confirmed by the 
same variation in /3„c with frequency: p-modes have f}„( close 
to 1 while mixed modes have f3„( close to 1/2. Mixed modes 
enter the frequency domain by the low part and one can see 
that p„c is closer to 1 at high frequency where the g nature of 
the mixed modes is less pronounced. Because of the weight- 
ing by the eigenfunction, the values of the splittings yield 
a mean rotation (£l)/2n that is no longer dominated by the 
outer layer contribution, since the mixed modes have ampli- 
tudes in the central region where the rotation rate is larger. 
As a result, one obtains a rotation rate that is higher than the 
surface value but still much lower than the rotation rate of 
the central regions. 
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- The red-giant branch: When the outer convective region pro- 
gresses inward, its uniform rotation extends deeper as well. 
The surface rotation rate decreases with time essentially due 
to an increase in the stellar radius. The core accelerates, 
and in the intermediate region, where the H-burning shell 
lies, a sharp rotation gradient develops (see iPalacios et al.l 
2006). The frequency spectrum of red giants is composed 
of g-dominated mixed modes and mixed modes with nearly 
equal p and g character. When the star evolves up the red- 
giant branch, the number of g-dominated mo des increases 
and la rgely outnumbers the p-g mixed modes (Dupr et et al.l 
2009). These modes mostly probe the core rotation (Beck 
et al. l2012h iDeheuvels et al.(l2012t iMosser et al.1 l2012ah . 
For model # 25 ( the last shown in Fig. |4j>, the core rotation 
amounts to 251/iHz and the expected excited frequencies are 
about 90 yuHz. This indicates a ratio £ = 5.6 (Eq.l25li. 
At such rotation rates, the first-order perturbation is not rele- 
vant to computing the rotational splittings. Nonperturbative 
methods must then be used. 

However, observations show that several RGB stars have 
core rotati on rates that are muc h lower than predicted by our 
models. (Moss er et al ] l2012ah . Indeed the observed values 
amount to a few hundred nHz, whereas here we find sev- 
eral dozen to a few hundred fjHz. This leads to the conclu- 
sion that the central regions of our models rotate too fast and 
that rotationally induced transport of the kind included in our 
models is not efficient enough to slow down the central rota- 
tion of red-giant stars. 



6. Seismic tests of transport of angular momentum: 
slowing down the rotating core of red-giant stars 

As in the solar case, one can wonder whether another mecha- 
nism, such as internal gravity waves or magnetic fields, operates 
more efficiently to shape the rotation profile, either at the red- 
giant phase or before. However, several issues must be discussed 
before we can reach such a conclusion. 



6. 1. Uncertainties on stellar modeling and transport of 
angular momentum 

A possible reason for this discrepancy is some inaccuracy in the 
current physical description of stellar models. Indeed several as- 
sumptions enter the rotationally induc e d tran sport as prescribed 
by IZahnl (1 19921) and iMaeder & Zahnl d!998l) . Besides, several 
other uncertainties on the input physics of stellar models can 
also affect the resulting rotation profile to some extent. In what 
follows, we discuss the sensitivity of the rotation profile of red- 
giant models to several uncertainties. 



6.1 .1 . Effect of convective overshoot 

A 1.3M star has a convective core during the MS. Overshoot 
shifts the tracks on the HR diagram, imitating tracks with higher 
masses. As a result, at a given radius in the RGB, the core of the 
model with overshoot has had less time to spin up, since it started 
its contraction later in the evolution. An overshoot of Q.lHp 
reduces the core rotation rate by about 32 % when = 3.73/? Q . 

On the other hand, we need to reduce the mass of the model 
with overshoot by 0.04M o to reproduce the RGB at the same 
location in the HR diagram. In that case, the model with over- 
shoot and M* = 1 .26M Q has a central rotation rate (again, when 



R+ - 3.73/? Q ) that is similar to the central rotation rate of the 
model without overshoot and M+ = 1.3M . 

6.1 .2. Effect of the initial rotation state 

It is well known that, after the early stages on the MS, the surface 
rotation rate does not depend on the initial conditions; it only 
depends on th e magnetic braking law (Eqs. l9land[T0l see, e. g. 
Maeder 2009 and references therein). We found that it is indeed 
the case for our models. The internal rotation profile down to the 
center is also independent of the initial conditions after the early 
stages of the MS. 

6.1 .3. Impact of the magnetic braking law 

We computed an evolution with the braking law of K88 with a 
value of Kw twice the solar calibrated value to reduce the rota- 
tion rate. We found that the surface rotation rate is slowed down, 
as expected, but the central rotation rate changes only slightly. 
For instance, at the TAMS, the surface rotation rate changes by 
20 %, and by 3 % at the base of the RGB, while the central ro- 
tation rate changes by 2 % at the TAMS, 3 % at the base of the 
RGB. 

We calculated another evolutionary sequence with the RM12 
magnetic braking law. With a coefficient Krm calibrated so as 
to yield the solar rotation rate in a solar model, our 1.3M is 
much more efficiently slowed down compared to the evolution 
computed with K88. The central rotation rate, however, changes 
by at most 2 %. 

6.1.4. Stability of the rotation profile 

The changes on the surface rotation rate caused by the mech- 
anisms described above do not change the central rotation rate 
significantly because the extraction of angular momentum from 
the core is not efficient enough. The core and the surface are 
not sufficiently coupled. One potential way of coupling the core 
to the surface is through instabilities that might arise because of 
rotation gradient becoming very steep (as seen in Fig. ITOb . The 
Rayleigh criterion requires that r 2 D. increases outwards. For our 
1 3M Q model, the rotation profile becomes unstable according to 
the Rayleigh criterion in two places, at the shell source and just 
below the convective envelope. However, the stabilizing effect of 
the density stratification overcomes this instability, as expressed 
by the Solberg-Hoiland criterion: 



N 2 = Nj + N 2 + > 0, 



(29) 



where is the Rayleigh frequency in a rotating medium, given 
at the equator by 



(30) 



According to the Solberg-Hoiland criterion, the rotation 
profiles we obtain are stable throughout the evolution. How- 
ever, t he Goldreich-Schubert-Frick e instability (GSF, see Fricke 



Goldreich & Schubert 1967) may occur in stellar condi- 
Hirschi & Maeder] d2010h have studied its effects on pre- 



1968; 
tions. _ 

supernova models and find that if < the GSF instability is 
always present regardless of the yu- and T-gradients. We imple- 
mented their prescription in CESTAM. Where Ni < 0, turbulent 
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Fig. 11. Top: Rotation profiles for 1.3M models at the end of the MS 
calculated assuming standard viscosity coefficients (standard), a verti- 
cal turbulent viscosity £) v computed with Ri r = 1 (Ri r ), a horizontal 
viscosity coefficient 100 times the standard value (A,), and all of the 
above (All). Bottom: The same for models at the base of the RGB when 
R t = 3.73^ (Av = 21.3/ii/z). 
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Fig. 12. Top: HR diagram showing the evolutionary track of a 1.2M G 
model computed with an overshoot of 0.1H P , a vertical turbulent vis- 
cosity Dv computed with Ri c = 1, a A, increased by a factor 10 2 (red 
line, dashed), and a 1.3M Q model computed without overshooting and 
with the "standard" viscosities (blue line, full). Dots indicate the loca- 
tion of models with Av = 23.3yuHz. Bottom: Rotation profiles for the 
models above. 



transport by the GSF instability operates with a viscosity coeffi- 
cient given by the solution of their Eq. (20), 

(N 2 T +N 2 +N 2 n )x 2 + 
[N 2 D h + N 2 (K + AO + N 2 n (K + 2D h )] x 

+N 2 n (D h K + D 2 ) = 0, (31) 

where K is the thermal diffusivity given by Eq. ( IA.7I ). The vis- 
cosity coefficient associated to the GSF instability is then given 
by Dgsf = 2x. 

We found that < in two regions from the end of the MS, 
just below the border of the outer CZ and at the shell source. At 
the shell source, the Q-gradient is high because of the contrac- 
tion of the layers left behind by the movement of the shell source, 
while the descending CZ lowers Q just above its border. In our 
models, the GSF turbulent viscosity is on the order of Dy in the 
unstable regions. The turbulent transport by the GSF instability 
reduces the Q-gradient in the unstable region at the shell source, 
reducing the central rotation rate by 3 % in models at the base of 
the RGB. The turbulent transport by the GSF instability does not 



have time to change the Q-gradient below the CZ before the 
convective instability sweeps over this region. 

6.1 .5. Uncertainties on the turbulent viscosity coefficients 

To extract more angular momentum, we can either increase the 
meridional circulation or increase the vertical turbulent viscosity. 
To increase the meridional circulation, one option is to decrease 
the inhibiting effect of A (see the expression for the yu-currents, 
Eq. IA.5l >. According to Eq. {§), an increase in will lead to a 
decrease in A. The prescription for the horizontal coefficient of 
turbu lence £>h i s still open to di scussi on. Prescriptions de rived 
by IZahnl (1 1992b . iMaederl (120031) . and iMathis et all d2004 can 
differ by as much as two orders of magnitude. We adopted, as 
a test, a value of that is 10 2 higher than that previously used 
(Eq.ES. 

The value of the critical Richardson number in Eq. ( 1B.4I I is 
usually assumed to lie between 1/4 and 1/6. Briiggen & Hille- 
brandt 120011) and lCanutol d 2002l) indicate, however, that we can 
expect instability for Ri 1. To increase the efficiency of the 
vertical turbulent viscosity in transporting angular momentum 
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Fig. 13. Top: Rotational splittings Sv„ e for £ = 1 modes as a function 
of the normalized frequency y/Av associated with the rotation profile 
shown in Fig. [T2]for the 1.2M G model. Middle: Corresponding Sv/B„ e 
for t = 1 modes as a function of the normalized frequency v/Av. Bot- 
tom: The same for £ = 2 modes. 



from the core, we computed an evolutionary sequence using 
Ri c = 1 to test the effect of the uncertainties. Figure QT| shows 
profiles of the rotation rate computed with different assumptions 
of the diffusion coefficients, at the TAMS and at the base of the 
RGB when = 3.73/? Q . Models have the same mass and ra- 
dius, hence the same Av. As expected, with an increased value 
of Ri c from 1/6 to 1, the central rotation rate decreases by 18% 
at the TAMS and by 30% at the base of the RGB. With a D h in- 
creased by two orders of magnitude (and Ri c = 1/6), the central 
rotation rate decreases by 84% at the TAMS and the base of the 
RGB. When both changes are included (an increase of Ri c and 
Dh as above), the central rotation rate is decreased by an order of 
magnitude, while we need a decrease of two orders of magnitude 
to reproduce the observed rotational splittings. 

6.1.6. A slowly rotating red-giant model 

We eventually computed an evolutionary sequence where we 
combined all effects causing a decrease in the rotation rate in the 
central region; an overshoot of Q.lHp and a vertical turbulent 
viscosity Dy computed with Ri c = 1, together with Dh increased 
by a factor 10 2 . This model was computed with M* = 1.2M 
because during the RGB the evolutionary track on the HRD of 
a l.2M computed with all these effects lies close to the evo- 
lutionary track of the 1.3M models used as comparison. The 
resulting rotation profile for models at the base of the RGB with 
same Av is significantly decreased in the central region with the 
central value £2 C = 21.5// Hz and is displayed in Fig. [12] 

Calculations of linear rotational splittings shows that the 
maximum splittings amounts to 5-6 ju Hz for model #20 at the 
base of the RGB (Fig. [T3~l >. This corresponds to a maximum 
value for 5v/j3„i of 12//Hz. The rotation gradient and the value 
of the central rotation are now too low for nonperturbative calcu- 
lations to give rise to significant corrections to the values of the 
linear rotational splittings. Such rotation splittings lead to a core 
rotation rate that is closer to the observations but still too high 



by almost one order o f magnitude compared with what is seen 
in recent observations dBeck et alJlioT 2; Deheuvels et a 
lMosseretaT]|2012ah . 

The seismology of red-giant stars emphasizes that an addi- 
tional mechanism is needed to achieve a slower core rotation of 
red-giant stars. 

7. Conclusions and perspectives 

We have computed stellar evolution models of low-mass stars, 
taking the transport of angular momentum by turbulent viscosity 
and merid ional circulati o n in r adiat ive zones into account. A s 
found by iPalacios et al] d2006l) and lEggenberger et alj d2010h . 
the core accelerates very rapidly in our models during the sub- 
giant and giant phases. We showed that the linear approximation 
for computing rotational splittings remains nevertheless valid 
for stochastically excited solar-type modes, for all evolutionary 
stages until the RGB. 

For red-giant stars, interpretations of recent Kepler observa- 
tions lead to a ratio between the core and surface rotation rates 
of five to ten, far lower than in our models where the ratio be- 
tween the core and surface rot ation rates approaches 10 3 . Similar 
conclusions were reached by Eggenberger et al.l d2012h . which 
was published during the submission process of the present pa- 
per. Nonperturbative calculations lead to complex frequency 
spectra, with in particular non symmetric multiplets ( see for in- 
stance without Coriolis force Lignieres et al. 2006; Ouazzani 
et al. 120 121 : 1 Rees e et al . 2006, with both centrifugal and Coriolis 
force). Observations show that several red-giant stars do indeed 
have such complex spectra. Their core rotation is fast enough to 
have entangled rotational splittings and mixed-mode spacings. 
In such cases, rotation may have to be studied with nonpertur- 
bative methods. On the other hand, other red-giant stars (Beck 
et al. 1201 21) show frequency spectra where the rotation splittings 
are easily identified as symmetric patterns around axisymmetric 
modes. The values of the corresponding splittings are quite low, 
close to or smaller than 1//Hz, and using the linear approxima- 
tion to derive the rotation splittings from stellar models is fully 
justified. 

We have computed the evolution of a stellar model including 
rotationally induced transport taking uncertainties in the values 
of the parameters entering such a description into account. We 
achieved a decrease in the central rotation rate of a red-giant star 
by an order of magnitude at most, mostly due to an increase in 
Dh by two orders of magnitude. The modification of the pa- 
rameters were chosen so as to enhance the transport of angular 
momentum so that the core rotation is slowed down as much 
as possible. The models in the red-giant phase were then slowed 
enough that the calculations of linear rotation splittings are valid. 
They were found to be one order of magnitude larger still than 
observed. This indicates that extraction of angular momentum 
from the core in our models is not efficient enough. 

We need to decrease it by another order of magnitude to re- 
produce the observations. A similar situation is encountered in 
the case of the Sun. Several possibilities have been proposed, 
the main candidates being magnetic fields and internal gravity 
waves. We showed that such mechanisms must operate or keep 
operating during the whole subgiant and red-giant phase. In this 
context, iTalon & Charbonnel d2008h show that transport of an- 
gular momentum by internal waves generated by the convective 
envelope can play an important role from the subgiant phase to 
the base of the RGB. They also find that this transport seems to 
have no major impact on stars ascending the RGB itself and later 
on. In the RGB phase, another mechanism must be called for to 
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explain the observations. This will be investigated further in fu- 
ture work. On the observational side, measurements of rotational 
splittings for subgiant stars are crucially needed to constrain the 
efficiency for their central layers to slow down. 

A standard, validated version of CESTAM will be available 
for download soon, together with grids of stellar models com- 
puted with several rotational velocities and adiabatic oscillation 
frequencies for selected models. 
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Appendix A: Meridional circulation 

The vertical component of the meridional circulation is given by 



s m — L/m the mean energy production inside a sphere of radius 
r, and p m = 3m/(4nr 3 ) the mean density. The thermal conduc- 
tivity x is given by 



4acT 3 



X = 



3i<p 

and the thermal diffusivity is 
X 4acT 3 



K = 



Hj 



pep 3Kp 2 Cp 
We also used 
dr 



dlnT 

the temperature scale height; 
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the logarithmic derivatives ofx', 
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w here c p , T, q, V, and V a d have t heir usual meanings (as defined 
in lKjp penhahn & Weigerill99lh . In Eq. (IA.lt above, M* is the 
reduced mass, 



! ) /Sine, 
st = I — — — I and s, 
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the logarithmic derivatives of e; and 
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(A. 10) 
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where m and p m are the mass and mean density inside an iso- 
bar, respectively, and Eg and E u denote the so-called £2- and 
/^-currents (following Palacios et al. 2003). The quantities ip and 
6 are obtained from the equa tion of state, and are defined as in 
iKippenhahn & Weigertl d 19911) . 



Appendix B: Turbulent coefficients of diffusivity 

The coefficient D er r is given bv (IChabover & Zahnll 19921) 



\d\aTj \<31nW/>,r 
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The Q- and ju-currents are given by 
Q? s + s g \g 
9 

P 



We use the prescription of M athis et alJ (120041) for the hori- 
zontal component of turbulent diffusivity D/,: 
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We take the coefficient p = 1 .5 x 10" 5 dRichard & Zahnll 19991) . 

For the vertical c omponent Dy, we use the prescription of 
lTalon&Zahr](fl997t) . 
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«p + ^stJ-Xv ~ -jjiXT + 1) 

In Eqs. (1A.4I > and dA.51 > we used the following quantities: s 
and Eg are the nuclear and gravitational energy generation rates, 



where Ri c = 1/6 is the critical Richardson number, K the ther- 
mal diffusivity, Nj and the chemical and thermal parts of the 
Brunt- Vaisala frequency, N 2 - N T + N 2 : 



jV2 = ^ (Vad _ V); 



(B.5) 
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